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Abstract. The VLT on Cerro Paranal was used to ob- 
serve four fields located at ~ 2.'3 from the center of the 
Galactic globular cluster NGC 6712 in the V and R bands. 
The resulting color-magnitude diagram shows a well de- 
fined main sequence reaching down to the 5er detection 
limit at V ~ 25, R ~ 23.5 or approximately 4 magnitudes 
below the main sequence turn-off, the deepest obtained 
so far on this cluster. This yields a main sequence lumi- 
nosity function that peaks at Mr ~ 4.5 and drops down 
to the 50% completeness limit at Mr ~ 8.5. Transfor- 
mation to a mass function via the latest mass-luminosity 
relation appropriate to this object indicates that the peak 
of the luminosity function corresponds to ~ 0.75 M©, a 
value significantly higher than the ~ 0.25 M©, measured 
for most other clusters observed so far. Since this object, in 
its Galactic orbit, penetrates very deeply into the Galactic 
bulge with perigalactic distance of ~ 0.3 kpc, this result 
is the first strong evidence that tidal forces have stripped 
this cluster of a substantial portion of its lower mass star 
population all the way down to its half-light radius and 
possibly beyond. 

Key words: globular clusters: general - globular clusters: 
individual: NGC 6712 - stars: luminosity function, mass 
function - stars: low-mass, brown dwarfs - stars: Popula- 
tion II 



1. Introduction 

NGC 6712 is a small (tidal radius ~ 8') and rela tively 
loose (c = 0.9) and faint (m-M = 15.6; Djorgovski |l993[ ) 



Galactic globular cluster (GC) that has not yet received 
much observational attention. Its main claim to fame so 
far is due to the presence in its core of the high luminosity 
X-ray burster X1850-086 whose optical counterpart may 
be a faint UV-excess object (Anderson et al. 1993). This 



fact presents somewhat of a puzzle since one would expect 
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such an X-ray source to be located in a highly concen- 
trated cluster where the stellar density favors its forma- 
tion via tidal capture of a neutron star (Hertz & Grind- 
lay 1985). Most other sources of this type have indeed 



been found in high density core collapse clusters suggest- 
ing that, perhaps, NGC 6712 has already undergone such 
an event in the past and is now in a state of re-expansion 



(Grindlay et al. |1988| ). 

This unusual situation may also be connected in some 
way to its Gala ctic orbit as computed recently by Daup- 
hole et al.(1996) that is fairly well restricted to the vicin- 
ity of the disk and penetrates very deeply in the Galactic 
bulge. This certainly means that one would expect this 
cluster to have undergone severe tidal shocking during the 
numerous encounters with both the disk and the bulge 
during its lifetime and the consequences on the dynami- 
cal status of the cluster to be significant and observable. 
A simple single-mass approximation of these effects was 



computed by Gncdin & Ostriker (1997) for both disk and 



bulge shocks under differing assumptions on the Galac- 
tic model with a resultant time to destruction as small as 
0.03 H . According to these calculations, then, the clus- 
ter should have evaporated long ago and at the very least 
may have lost a very substantial portion of its original 
mass during its lifetime. 

Clearly, this catastrophe should be well impressed on 
its present day distribution of stars on the main sequence 
(MS) with its lowest mass members beyond the half-mass 
radius particularly vulnerable to escape. This effect may 
well have been detected already in M4, another cluster 



at significant risk of tidal disruption (Kanatas et al. 1994; 
Pulone et al. 1998a), but until this cluster's structural pa- 



rameters are pinned down more accurately this remains 
still speculative. There is, therefore, much interest today 
in determining accurately the present day mass function 
(PDMF) of NGC 6712 to look for the signature of such 
powerful effects. Currently available observations of the 
color-magnitude diagram (CMD) of this cluster, however, 



only reach to just above the MS turn-off (Cudworth 1988 



Anderson et al. 1993) and are, therefore, of limited use for 



this task. In order to push the observations well into the 
relevant part of the MS below the turn-off (TO), the VLT 
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was used to probe deeply into this cluster with its unprece- 
dented sensitivity and resolution. This paper describes the 
first results of these observations that give clear evidence 
that there is indeed a distortion of the MF of NGC 6712 
with respect to that of other dynamically much less dis- 
turbed clusters. 



2. Observations and Data Analysis 

The observational data used in this paper were col- 
lected duri ng th e Science Verification (SV; Leibundgut 
& Renzini 1998) phase of the first 8m-diameter Very 
Large Telescope (VLT) at ESO, using the VLT Test 
Camera (VLT-TC). Readers interested in the VLT and 
its instruments should consult ESO's world-wide web at 

while the 



tittp : //www . eso . org/paranal 



scope of the 

VLT Science Verification is described in Leibundgut, 
De Marchi, & Renzini (199S). Images of the globular clus- 



ter NGC 6712 were taken with the VLT-TC in the V and 
R bands. With a 2, 080 x 2, 048 square pixel detector and 
a plate scale of 0'.'045 pixel" 1 , the VLT-TC offers a field of 
view of ~ 90" x 90". SV observations, however, were ob- 
tained with an electronically enforced 2x2 binning of the 
CCD, so that the actual size of each pixel in these images 
is of 0"091 on a side. Observations of four regions of the 
cluster are available, located between one and t wo times 
the half-light radius (r M ~ 78"; D,jorgovski |l993] ). The co- 
ordinates (J2000) of the center of each field are given in 
Table 1 along with the total exposure time in each band. 
Fields F 1 and F2 were observed during the night of 1998 
Aug 23, and F3 and F4 on 1998 Aug 27. 

Table 1. VLT-SV observations of NGC 6712 
Field RA DEC r t v t R 



Fl 
F2 
F3 
F4 



18:53:09.1 
18:53:02.0 
18:53:10.0 
18:52:59.5 



-8:41:17 
-8:41:02 
-8:43:04 
-8:43:25 



145" 
105" 
130" 
135" 



900 s 1800 s 
1800 s 
2700s 2467s 
900 s 900 s 



In our investigation we have used the standard SV 
combined datasets corresponding to Field 1, 3, and 4 which 
are shown in Figure 1 (Field 2 was not used as V-band 
observations are not available). The quality of these im- 
ages is excellent, with seeing full width at half maximum 
(FWHM) always of order ~ 0'.'6. As can be seen in Fig- 
ure 1, however, the stellar density increases considerably 
towards the cluster center, thus making it progressively 
more difficult to accurately measure the magnitude of faint 
stars. Since the robustness of the luminosity function (LF) 
that one could determine from such images is inversely cor- 
related with the level of crowding, we have restricted our 
analysis to the quadrants farthest away from the center 
(see Figure 1), where the star density and the number of 
bright objects is smaller. The distance between the center 
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Fig. 1. Negative R-band image of fields F 1, F3, and F4 
(top to bottom). North is up and East points to the left. 
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of each quadrant and the nominal cluster center is given 
in Table 1. 

The IRAF automated star detection routine ap- 
phot.daofind was applied to the data, by setting the de- 
tection threshold at 5 a above the local background level. 
We then carefully examined by eye each individual object 
detected by daofind and discarded heavily saturated stars 
and a few extended objects whose FWHM exceeds by a 
factor of two or more that typical of stellar sources in our 
frames (FWHM ~ 0'.'6). We identified in this way 328, 
500, and 568 well defined stellar objects in both bands 
respectively in the outer quadrants of Field 1, 3, and 4, 
and measured their fluxes using the standard IRAF ap- 
phot.phot aperture photometry routine, by setting the ob- 
ject radius to r — 0'.'46 and the background annulus from 
0'.'46 to 0'.'73. We have estimated that this choice of aper- 
ture and background annulus samples a fraction of the 
total stellar flux (defined as that falling within an aper- 
ture of 6") which, on average, corresponds to er ~ 0.48 
and £r ~ 0.45 respectively in the V and R bands. In- 
strumental magnitudes were converted into the standard 
Johnson system using the calibration provided by ESO as 
part of the SV data release, namely: 

mag = —2.5 log (cs/t) — C\ — C2 colorterm — C3 airmass(l) 

where c is the number of counts measured within the se- 
lected aperture, e is the corresponding encircled energy, 
and t is the exposure time. The coefficients Ci, C2, and C3 
are given in Table 2. 

Table 2. Photometric coefficients 



Night 


Band 


Cl 


C2 


C3 


23 Aug 


V 


-26.764 


0.057 


0.187 


23 Aug 


R 


-26.916 


0.054 


0.122 


27 Aug 


V 


-26.798 


0.064 


0.210 


27 Aug 


R 


-26.980 


0.059 


0.168 



The resulting CMD is shown in Figure 2, where the 
data from Field 1, 3, and 4 are marked with boxes, trian- 
gles, and crosses, respectively. Stars brighter than R ~ 17 
are heavily saturated and are not plotted. Objects brighter 
than R ~ 18 are likely to be affected by some degree of non 
linearity in the detector response and, as such, their mag- 
nitudes are not reliable, while stars fainter than R ~ 18.5 
are comfortably within the linear regime of the camera. 

We have compared the V and R magnitudes of sev- 
eral stars in Field 1 and 3 with those measured on a set of 
short WFPC 2 exposures in F555W and F675W extracted 
from the HST archive. We find that the zero point of our 
photometry agrees with that of the HST (VEGAMAG sys- 
tem) to within ~ 0.1 and ^0.2 magnitudes respectively 
in the V (F555W) and R (F675W) bands. The still pre- 
liminary calibration of the VLT-TC coupled with the un- 
certainty in our aperture corrections because of crowding 
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Fig. 2. Color-magnitude diagram of the stars detected in 
Fields Fl (boxes), F3 (triangles), and F4 (crosses). The 
solid line shows the location of the MS as predicted using 
the models of Baraffe et al. (see text) 



prevent us from determining the zero point of our pho- 
tometry with greater accuracy. Nevertheless, because the 
VEGAMAG system does not reflect exactly the Johnson 
system (particularly F675W), we consider this agreement 
very good. We should also point out that the results of 
the investigation presented in this paper are insensitive to 
errors of a few tenths of a magnitude in the zero point (see 
Section 4 below). 

The MS of NGC 6712 is rather well defined from the 
turn-off (TO) at R ~ 19 down to R ~ 23, where it 
broadens due to the increasing photometric errors. In fact, 
the accuracy of our measurements varies from less than 
- 0.05 mag at R ~ 18.5 to - 0.35 mag at R ~ 23. The 
solid line plotted over the CMD in Figure 2 corresponds 
to the expected location of the MS a s pred icted using the 
theoretical models of Baraffe et al. ( 1997 ) for the metal- 
licity, distance, and reddening a ppropr iate to NGC 6712 
([Fe/H]= -1.0, from Zinn & West |l984j ( m-M ) = 14.16 



and E(B - V) = 0.46, from Djorgovski |l993j ). The good 
agreement between our data and the models suggests that 
the latter are reliable, and that the use of the correspond- 
ing M-L relations when converting a MF into a LF (see 
Section 4) should give accurate results. The CMD in Fig- 
ure 2 also reveals a few objects on both sides of the MS 
which are likely contaminating field stars due to the low 
galactic latitude of the cluster (bn ~ —4.3 at ~ 500 pc 
below the Galactic plane). 
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3. The Luminosity Function 

From the CMD of Figure 2, we have measured the LF of 
MS stars in each field by counting the number of objects 
in 0.5 mag bins along the R axis. In our exercise, we have 
not accounted for the contamination due to field stars and, 
as a result, our LFs are probably an overestimate of the 
true distribution. This effect, however, becomes significant 
only below R ~ 21, where photometric incompleteness 
increases considerably and is likely to represent the largest 
source of uncertainty in our measurements. In fact, the 
artificial star tests that we have carried out to estimate 
the level of completeness as a function of magnitude (see 
Table 3) show that at R ~ 22 we only sample two thirds of 
the total population, and that this fraction drops to 50 % 
at R ~ 23.5. 



Table 3. Photometric completeness as a function of mag- 
nitude 



Rmag — 


18.0 


19.0 


20.0 


21.0 


22.0 


23.0 


23.5 


Fl 


1.00 


0.92 


0.83 


0.77 


0.70 


0.55 


0.48 


F3 


1.00 


0.98 


0.90 


0.82 


0.60 


0.50 


0.40 


F4 


1.00 


0.95 


0.88 


0.78 


0.70 


0.62 


0.58 



The LFs measured in this way and corrected for photo- 
metric incompleteness are shown in Figure 3 as a function 
of the R-band magnitude. The LFs have been registered 
through a vertical shift in the logarithmic plane by impos- 
ing a least square fit in the range 19 < R < 23. The error 
bars associated with each point reflect the poisson statis- 
tics of the counting process (and include the correction for 
incompleteness). These LFs can be directly compared to 
one another as they have all been measured at the same 
radial distance from the center (r ~ 2.'3 or ~ 1.7 times the 
half-light radius rhi). They show the same overall trend, 
i.e. an increase with decreasing luminosity up to a peak at 
R ~ 20 (close to the TO luminosity), and from there they 
all flatten out and possibly drop with decreasing luminos- 
ity even after the incompleteness of our photometry has 
been accounted for. And indeed, to ensure that our LFs are 
robust, we have not included in Figure 3 any datapoints 
whose associated photometric completeness is worse than 
50 %. 



4. Discussion and Conclusions 

Stars brighter than R ~ 19 have already evolved off the 
MS and, therefore, their LF provides no information on 
the underlying MF without uncertain corrections for evo- 
lution (Scalo |l99"8 ). Moreover, because of saturation at the 
bright end of our CMDs, the brightest portion of our LFs 
is uncertain. For cluster stars which are still on their MS, 
however, the LFs in Figure 3 directly reflect the PDMF of 
the local population and immediately indicate a relative 
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Fig. 3. Luminosity functions of MS stars measured in 
Fields Fl (boxes), F3 (triangles), and F4 (crosses), af- 
ter correction for photometric incompleteness 



deficiency of low mass objects with respect to the stars 
with the TO mass (~ 0.8 M ,), as we discuss below. 

Indeed, the most important conclusion that one can 
draw from Figure 3 is that the shape of the LFs completely 
deviates from that of any other GC for which relatively 
deep photometric data are available near the half-mass ra- 
dius. Observations carried out with the WFPC 2 on board 
the HST o ver th e past few years (Pare sce, D e Marchi & 
Rom anicllo |l995| ; Cool, Piott o, fc K i ng |l996| ; Elson et al. 



1995 



al 
that 



1995b 



1996a, 1997 



De Marchi & Paresc e |l995a 
Pio tto, Cool, & Kin g |1997| ; Pulone et al. [19984 Kin S et 
1998; De Marchi 1998) have consistently revealed LFs 
near the cluster half-mass radius, increase with de- 
creasing luminosity from the TO magnitude all the way 
down to about My — 11 (— 0.25 M Q ,) where they flatten 
out and drop at fainter luminosities. Inverted LFs such 
as those shown in Figure 3 have been observed right in 
the core of high density GCs (47Tuc, NGC6397, M 15) 
but in those cases a simple isothermal model of a cluster 
in equilibrium can easily explain this effect as being due 
to mass segregation (Paresce, De Marchi, & Jedrzejew- 
ski 



1995 



King, Sosin, & Cool 1995; De Marchi & Paresce 
1996b). More complete multi-mass King-Michie models 
show, however, that thermal relaxation is much less effi- 
cient (if at all) at depleting low-mass stars near the half- 



mass radius (see Pulone, De Marchi, & Paresce 1998b), 
and we cannot therefore trace the origin of the LFs that 
we observe back to the effects of mass segregation alone. 
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To make it easier to compare the LF of NGC6712 
with that of other clusters, we display it in Figure 4 as 
a function of the absolute R-band magnitude, assuming 
(m - M) = 14 16 and E(B - V) = 0.46 or A R = 1.15 
(Djorgovski 1993). Rather than showing the three individ- 
ual LFs, we have combined them together into one single 
function by averaging their values in each magnitude bin, 
and have taken the standard deviation as a measure of the 
associated uncertainty (error bars) . The dashed line shown 
in Figure 4 corresponds to the LF of the low-metallicity 
cluster NGC 6397 as measured by King et al. (1998), while 
the dot-dashed line reproduces th e LF of the metal rich 
cluster 47Tuc from Hesser et al. ( |l987| ). Both LFs have 
been translated into the R-ba nd by using the M-L rela- 
tionship of Baraffe et al. ( 1997 ) for the appropriate metal- 
licity, i.e. the magnitude corresponding to each observed 
point in the LF has been converted into a mass which has 
then been used to read the corresponding magnitude in 
the R band from the appropriate M-L relation. The size 
of each magnitude bin has also been rescaled to reflect 
the difference in the slopes of the M-L relationships for 
different bands. We have selected NGC 6397 and 47 Tuc 
as they both have accurate LF measurements at and be- 
low the TO luminosity, where we have normalized them to 
our observations, and because the metal content of these 
clusters nicel y bra ckets that of NGC 6712 QFe/H]= -1.0; 
Zinn & West 1984). It should, nevertheless, be clear that, 
due the uncertainties in the theoretical M-L relations and 
in the observed LFs, our comparison will only provide an 
indication of the true differences. 

The difference between these two LFs and that of 
NGC 6712 is striking. While the LFs of NGC 6397, mea- 
sured at 1.6 rm, shows a steep increase starting from the 
TO, the LF of NGC 6712 sampled at ~ 1.7r hl slowly 
drops from the TO all the way to the detection limit at 
Mr — 8.5. We would like to point out that the discrep- 
ancy is so large that to bring the two LFs into agreement 
would require us to have underestimated the photometric 
incompleteness by a factor of ~ 10. The same reasoning 
holds true for the LF of 47 Tuc, which has been measured 
at ~ 5rhi- This difference must thus be physical and re- 
flect the properties of the local stellar population. 

Under the simple assumption that the MF should be 
represented by an exponential distribution in the mass 
range 0.4 — 0.8 M©, (a reasonable hypothesis given the 
narrow mass range), we have used the M-L relationship 



2.5 



of Baraffe et al. ( 1997 ) appropriate for the metallicity of 
NGC 6712 to reproduce the observed LF. We obtain a 
fairly reasonable fit to the observations with a power-law 
distribution of the type dN/dm oc m 15 (Salpeter's IMF 
would be dN/dm oc m -2 35 ), in which the number of ob- 
jects decreases with mass (solid line in Figure 4). 

Richer et al. ( |1991| ) and, more recently, De Marchi & 
Par esce (1 997), Vesperini & Heggie (1997), and Pulone et 
al. ( |l998b!) have convincingly shown that near the cluster 
half-mass radius the LF should closely reflect the IMF, 
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Fig. 4. Boxes: average of the three LFs shown in Figure 3, 
converted to absolute magnitude (Mr); dashed and dot- 
dashed lines: LFs of NGC 6397 and 47 Tuc translated to 
the R band and normalized to ours at the TO (Mr ~ 4); 
solid line: best fitting power-law MF (dN/dm oc m ) 



as dynamical modifications should leave these regions al- 
most untouched. In fact, the internal relaxation mecha- 
nism governed by energy equipartition through two- and 
three-body encounters mostly affects the region within 
a few core radii, while the interaction with the Galactic 
tidal field is expected to simply speed up the evaporation 
of light stars near the tidal boundary, but none of these 
processes should, in principle, significantly alter the prop- 
erties of stars located close to the much safer half-mass 
radius area. 

If this were true for NGC 6712 as well, one should con- 
clude that this cluster is the only one so far to feature 
an inverse IMF (increasing with mass) that has not been 
observed in any other environment. While this hypothesis 
cannot be safely ruled out, there are far better reasons 
to believe that NGC 6712 might have experienced a much 
stronger interaction with the Galaxy than any other of 
the clusters studied so far. And indeed, with a perigalac- 
tic distance smaller than 300 pc this cluster ventures so 
frequently and so deeply into the Galactic bulge (Daup- 
hole et al. 1996) that it is likely to have undergone se- 
vere tidal shocking during the numerous encounters with 
both the disk and the bulge during its lifetime. The lat- 
est Galactic plane crossing could have happened as re- 
cently as 4 10 6 year ago (Cudworth 1988| ) which is much 
smaller than its half-mass relaxation time (5 10 8 yr). Such 
an event might have imparted strong modifications on the 
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mass distribution not only of the stars in the cluster pe- 
riphery but also well into its innermost regions, perhaps 
even reaching the core where it could have triggered a pre- 
mature collapse because of tidally induced relaxation (see 
Kundic & Ostriker H997I and Gnedin & Ostriker |1997| for 



a detailed description of this mechanism). 

As a result of such a catastrophe, it would be surpris- 
ing if the present-day MF were still to bear any mem- 
ory of its parent IMF anywhere in the cluster, including 
the half-mass radius region. Vesperini & Heggie (1997) 
have estimated that these effects would substantially de- 
crease the slope of a simple power law MF, much in the 
same way as we are observing here. We, therefore, con- 
clude that the VLT has revealed the consequences of the 
strong tidal stripping that the Galaxy (and particularly 
its bulge) exerts on GCs orbiting close to the center, and 
which might have contributed to the destruction of an ini- 
tially much more numerous population of GCs (Aguilar, 
Hut, fc Ostriker |1988|; Vesperini |1997[). Although Kanatas 



et al. (1994) and Piotto et al. (1997) had speculated that 



similar events could have happened respectively in M4 
and NGC6397, the result that we show here is the first, 
clear, unambiguous detection of this mechanism. To char- 
acterize the strength and extent of these phenomena more 
accurately would require the investigation of the MS popu- 
lation outside the half-mass radius in many more clusters, 
and possibly at larger distance from the Galactic center, 
so as to probe the intensity of the stripping process as a 
function of the depth of the Galactic potential well. If the 
Z component of the space velocity of NGC 6712 is indeed 
appropriate for a halo cluster, as suggested by Cudworth 
( |1988| ), then this violent stripping process might not be 
restricted only to objects orbiting the innermost Galactic 
regions. 

Acknowledgements. We are indebted to the referee for critical 
comments that have greatly improved our paper. 
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